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Abstract. From the spectra and light curves it is clear that SNIa events are 
thermonuclear explosions of white dwarfs. However, details of the explosion 
are highly under debate. Here, we present detailed models which are con- 
sistent with respect to the explosion mechanism, the optical and infrared 
light curves (LC), and the spectral evolution. This leaves the description of 
the burning front and the structure of the white dwarf as the only free pa- 
rameters. The explosions are calculated using one-dimensional Lagrangian 
codes including nuclear networks. Subsequently, optical and IR-LCs are 
constructed. Detailed NLTE-spectra are computed for several instants of 
time using the density, chemical and luminosity structure resulting from 
the LCs. The general methods and critical tests are presented (sect. 2). 

Different models for the thermonuclear explosion are discussed including 
detonations, deflagrations, delayed detonations, pulsating delayed detona- 
tions (FDD) and helium detonations (sect. 3). Comparisons between theo- 
retical and observed LCs and spectra provide an insight into details of the 
explosion and nature of the progenitor stars (sect. 4 & 5). We try to answer 
several related questions. Are subluminous SNe la a group different from 
'normal' SN la (sect. 5)? Can we understand observed properties of the LCs 
and spectra (sect. 4)? What do we learn about the progenitor evolution and 
its metallicity (sect. 3, Figs. 4,5)? Do successful SN la models depend on 
the type of the host galaxy (Table 2)? 
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Using both the spectral and LC information, theoretical models allow for 
a determination of the Hubble constant independent from 'local' distance 
indicators such as 5 Cephei stars. is found to be 67 it 9 km s~^Mpc~^ 
and, from SN1988U, Qo equals 0.7 ± 1. within 95 % confidence levels. 



1. Introduction 

During the last few years it became evident that Type la supernovae are a 
less homogeneous class than previously believed (e.g. Barbon et al., 1990, 
Pskovskii 1970, Phillips et al. 1987). In particular, the observation of sub- 
luminous SN1991bg (Fillipenko et al. 1992, Leibundgut et al. 1993) raised 
questions on the SN-rate and whether we have missed a huge subgroup. 
The possible impact on our understanding of supernovae statistics and, 
consequently, the chemical evolution of galaxies must be noted. 

It is generally accepted that SNe la are thermonuclear explosions of 
carbon-oxygen white dwarfs (WD) (Hoyle &; Fowler 1960). However, details 
of the scenario are still under debate. For discussions of various theoretical 
aspects see e.g. Wheeler & Harkness (1990), Woosley Sz Weaver (1994), 
Canal (1995), Nomoto et al. (1995), Wheeler et al. (1995). 

What wc observe as a supernovae event is not the explosion itself but 
the light emitted from a rapidly expanding envelope produced by the stellar 
explosion. As the photosphere recedes, deeper layers of the ejecta become 
visible. A detailed analysis of the LCs and spectra gives us the opportunity 
to determine the density, velocity and composition structure of the ejecta. A 
successful application of these constraints, however, requires both accurate 
early LC and spectral observations and detailed theoretical models. While 
SN la LC data have enormously increased, until recently (Harkness 1991; 
Hoflich et al. 1991), models are often hampered by inadequate physical 
assumptions like a constant opacity, crude gamma-ray deposition schemes 
and a simplified treatment of the ionization balance, neglect of line scat- 
tering and by the use of the diffusion approximation (e.g. Livnc & Arnctt 
1995). Based on our detailed models which overcome the approximations 
just mentioned, we have investigated the validity and influence of the phys- 
ical assumptions made in LC and spectral calculations and estimated the 
accuracy of our models. The importance of a consistent treatment of the ex- 
plosion mechanism, LCs and spectra became evident. Some of the tests and 
the basic outline of our approach is given in section 2. In section 3, the basic 
explosion mechanism are discussed. After investigating the general proper- 
ties of our LCs and spectra, individual comparisons with observations are 
presented. Finally, we address the question on Hg and Qo- 
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2. Brief Description of the Numerical Methods 

2.1. HYDRODYNAMICS 

The explosions are calculated using one-dimensional Lagrangian hydro with 
artificial viscosity (Khokhlov, 1991ab) and radiation-hydro codes including 
nuclear networks (Hoflich & Khokhlov 1995). The latter code is based on 
our LC code that solves the hydrodynamical equations explicitly by the 
piecewise parabolic method (Collela and Woodward 1984) and includes the 
solution of the radiation transport implicitly via the moment equations, 
expansion opacities, and a detailed equation of state. Typically, 300 to 500 
depth points are used. Radiation transport has been included to provide 
a smoother transition from the hydrodynamical explosion to the phase of 
free expansion. We omit 7-ray transport during the hydrodynamical phase 
because of the high optical depth of the inner layers. Nuclear burning is 
taken into account using Thielemann's network (Thielemann, Nomoto & 
Hashimoto 1994 and references therein). During the hydro, an a-network 
of 20 isotopes is considered to properly describe the energy release. Based 
on a network of 216 isotopes, the final chemical structure is calculated by 
postprocessing the hydrodynamical model. The accuracy of the energy re- 
lease has been found to be about 1 to 3 % in the reduced network compared 
to about 15 % for the equation net work used previously (Khokhlov, 1991) 
which, in general, overestimate the expansion velocity accordingly. 

2.2. LIGHT CURVES 

Based on the explosion models, further hydrodynamical evolution, and 
bolometric and monochromatic light curves are calculated using a scheme 
recently developed and widely applied to SN la (Hoflich et al. 1995 and ref- 
erences therein). In principle, it is the same code as that described above, 
but nuclear burning is neglected and 7 ray transport is included via a Monte 
Carlo scheme. In order to allow for a more consistent treatment of scat- 
tering, we solve both the (two lowest) time-dependent radiation moment 
equations for the radiation energy and the radiation flux, and a total energy 
equation (Mihalas 1978). The Eddington factors are obtained by formal in- 
tegration of the radiation transport equation. The opacities have been cal- 
culated under the assumption of local tlicrniodynamical equilibrium. This 
is a reasonable approximation since diffusion time scales are governed by 
layers of large optical depth. The Planck and the Rossland mean of the 
absorption opacity and extinction, respectively, have been calculated based 
on an approach similar to Karp et al. (1977) but transformed to the comov- 
ing frame. Relativistic terms are neglected. As in any Sobolev approach, it 
is assumed that the resonance region of lines is small. Therefore, the mean 
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Figure 1. Bolometric LC of the delayed detonation model N21 for various assump- 
tions (N21: all effects are included as described in the text; N21K2, N21K1 & N21K05: 
= Kp = const with 0.20, 0.1 and 0.05 cm^ /g, respectively; N21S: kp = Kp, N21A; 
extinction = true absorption; N21D: diffusion approximation). In addition the energy 
release due to the radioactive decay of ^'^'Ni and ^^Co is shown (dotted line). If constant 
opacities are assumed, explosion models cannot be tested because the strong temperature 
dependence of the opacity for T < 20000i^ is neglected. This effect determines the shape 
of LCs near maximum light. For reasonable variations in k (Shigeyama et al. 1993), the 
time of maximum light is shifted by up to 5 days and the maximum brightness may be 
off by about 50% . The shape of the LCs differ strongly. If the temperature dependence 
of the opacity is included, the spread in the results shrinks to about 20 % both with re- 
spect to the rise time and absolute flux. The slope becomes comparable within 10% if the 
radiation transport equation is solved rather than assuming the diffusion approximation. 




Figure 2. B and V LCs of N21 for various assumptions (see Fig. 1). Apparently, both the 
diffusion approximation (N21D) and the neglection of line scattering (extinction = true 
absorption, N21A) are inadequate. The the local cooling of the radiaton field is strongly 
overestimated since, in SNela, scattering exceeds the true absorption by a factor of 100 
to 10000 (Hoflich et al. 1993, Hoflich 1995). In comparison, the corresponding error by 
setting Kp * a = Hp in the energy equations results in errors of about 0.1 m because 
the thermal coupling differs by a much smaller factor (~ 5 ... 10). The surprisingly small 
sensitivity of the resulting LCs can be understood in terms of the energy conservation. 
Namely, during early times, the energy balance is determined by radioactive decay and 
the adiabatic expansion if the opacity exceeds a certain limit. 



free path of a photon along a certain ray is given by the product of inverse 
absorption probabihties weighted with the absorption probabilty of transi- 
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tions being closer to the comoving frame wavelength. The expansion effect 
enters both the Planck and the Rosseland mean but the former in second 
order, only, because the Planck mean is a linear average of k,^. The main 
problem with the Planck opacity is related to line scattering. In SNela, scat- 
tering exceeds true absorption by several orders of magnitude. Therefore, 
scattering must be included. Since 1993, we have used an equivalent-two 
level approach calibrated by detailed NLTE-calculations (see below, Hdflich 
1995). We found that the thermalization process is governed by redistribu- 
tion (i.e. collisions) of energy within the atomic fine structure levels and 
not by the scattering in about a dozen UV lines (Pinto, talk at this meet- 
ing). The former effect is dominant since coUisional cross sections depend 
exponentially on the energy difference between levels. In fact, the two-level 
approach underestimates the size of the thermalization by several orders 
of magnitude as it is well known in the field of stellar atmospheres (e.g 
Mihalas 1978, Ayres 1989). For details, see Hoflich (1995). 

Monochromatic colors B, V, R, and I are calculated using 100 discrete 
wavelength bands and formal integration along rays. For several models, 
detailed NLTE-spectra have been constructed (Hoflich, 1995). The colors 
based on the NLTE atmospheres and LC calculations have been compared. 
Typically, the disagreements between these two techniques in V, B, R, and 
I are below 0.05 to 0.1™ near maximum light and remain between 0.2 to 
0.4™ at late time. The error is small because most lines are scattering dom- 
inated. Consequently, the flux is redistributed within a wavelength range 
corresponding to the Doppler shift produced by the photospheric expan- 
son. However, this frequency shift is small compared to the bandwidth of 
the UBVRI-filter. With increasing time, the forbidden lines dominate the 
flux and, consequently, the error in our monochromatic colors increases. 
Another uncertainty when comparing with observations is due to the filter 
functions. Here, we use transmission functions given by Bessell (1990). 

In light of the uncertainties, we have compared our opacities with those 
based on approximations by Jeffery (1995). We found good agreement. 
Moreover, the sensitivity of the resulting LCs variations in the uncertain 
quantities have been studied in great detail (Hoflich et al. 1993, Khokhlov 
et al. 1994). Some of the results are shown in Figs. 1 and 2. Other, indepen- 
dent tests are provided by the comparisons with observations. We find that 
our bolometric correction is consistent with the observation of SN1992A, 
and the distances predicted by our models agree with those based on S- 
Ceph stars in IC4182, NGC5253 and 4536 (Hoflich k Khokhlov 1995). The 
former test probes the overall energy distribution, whereas the latter test 
is sensitive to errors in the absolute, monochromatic fluxes. 
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Figure 3. Various quantities for M36 15 days after the explosion. Distances are given 
as a function of wavelength at which the monochromatic optical depth reaches 0.1 and 
1, respectively (upper left), Temperature T, density p, energy deposition due to radioac- 
tive decay and Rosseland optical depth are given as a function of distance and, for 
comparison, Tgrey for the grey extended atmosphere (upper right), and chemical profile 
arc presented for the most abundant elements (lower left). No well defined photosphere 
exists. The need for spectral analyses consistent with respect to the explosion mechanism, 
7— ray transport and LC calculation is obvious. Already at maximum light, the chemi- 
cal distribution of elements does not follow the density profile but shows large individual 
variations in the line forming region, i.e. between 1 and 2 10^^ cm. The smallest expansion 
velocity indicated by line shifts are predominantly determined by the chemical profiles 
and not by the density. The energy deposition and excitation due to 7-rays becomes 
important within the photosphere. Consequently, the luminosity cannot be assumed to 
be depth independent for the construction of synthetic spectra. 

2.3. SYNTHETIC SPECTRA 

A modified version of our code for Nlte Extended ATmospheres (NEAT) is 
used. For details see Hofiich (1990, 1995) and references therein. 

Stationarity is assumed where the density, chemical profiles and the 
luminosity as a function of the radial distance r arc given by results of the 
hydrodynamical explosion and of the LC calculations including the Monte- 
Carlo scheme for the 7-ray transport (Fig. 3). 

The radiation transport equation is solved in the comoving frame in- 
cluding relativistic terms. For the determination of the first moment of the 
intensity, we solve the radiation transport equation for strong lines includ- 
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ing 'quasi' continua (weak lines are treated in a Sololev approximation, see 
below) in the comoving frame, i.e. 

OpI = x{S - 1) 

I = I{f^^r,uy, S = S{r,ij); x = x{r,i^) 

{li: cosine between the radial direction and I; r : radial distance; v(r) : 
radial velocity; S: source function; absorption coefficient; v. frequency). 
The operator Op can be written as follows (Castor, 1974): 



Op = 
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The term in front of ^ can be interpreted as the classical Doppler shift. 
The second terms of the partial derivatives with respect to r, jj, and the 
logarithmic derivative of v(r) correspond to the advection and aberration 
effects. For practical purpose, we solve the non rclativistic transport equa- 
tion using a Rybicki scheme (Mihalas et al., 1975, 1976). Overlapping lines 
in the comoving frame cannot be treated because the system is solved as 
a boundary problem in the frequency space which, from the concept, dis- 
allows a propagation of information toward higher /lower frequencies for 
expanding/collapsing envelopes. The same limitation should also apply to 
Nugent et al. (1995a) who use the same method. However, the problem of 
overlapping lines can be and has been solved (Hoflich 1990) in full anal- 
ogy with the problem of partial redistribution using a Fourtrier-like scheme 
(Mihalas et al, 1976b). Although our code can deal with partial redistri- 
bution, this effect is neglected to save CPU-time. A technical problem is 
the fact that each line must be sampled by at least 10 frequency points to 
provide a sufficient accurate representation of the S/Sv-term (Mihalas et 
al. 1975). Therefore, we must limit this approach to 2000 lines. 

Blocking by lines other than the strong ones is included in a 'quasi' 
continuum approximation, i.e. the frequency derivative terms in the radia- 
tion transport equation is included in the narrow line limit to calculate the 
probability for photons to pass a radial sub-shell along a given direction fi 
(Castor, 1974, Abbot and Lucy 1985, Hoflich, 1990). Note that the infor- 
mation on the exact location of the interaction within a subshell gets lost 
translating into an uncertainty in the wavelength location of ?a 10~^. 
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The statistical equations are solved consistently with the radiation trans- 
port equation to determine the non-LTE occupation numbers using both 
an accelerated lambda iteration (sec Olson ct al. 1986) and an equivalent- 
two level approach for transitions from the ground state which provides an 
efficient way to take the non-thermal fraction of the source function into 
account during the radiation transport and, effectively, accelerate both the 
convergence rate and the stability of the system. A comparison of the ex- 
plicit with the implicit source functions provides a sensitive tool to test for 
convergence of the system of rate and radiation transport equations. 

Excitation by gamma rays is included. Detailed atomic models are used 
for up to three most abundant ionization stages of several elements, i.e. 
(He), C, O, Ne, Na, Mg, Si, S, Ca, Fe taking into account 20 to 30 levels and 
80 - 150 transitions in the main ionization stage. Here, we use detailed term- 
schemes for C II, O II, Ne I, Na I, Mg II, Si II, S II, Ca II and Fc II or Fe III. 
The corresponding lower and upper ions are represented by the ground 
states. The energy levels and cross sections of bound-bound transitions are 
taken from Kurucz (1993). In addition to our NLTE-transitions, a total 
of « 300,000 lines out of a list of 31,000,000 (Kurucz, 1993) are included 
for the radiation transport. For these lines, we assume LTE population 
numbers inside each ion. To calculate the ionization balance, excitation by 
hard radiation is taken into account. LTE-line scattering is included using 
an equivalent-two-level approach calibrated by our NLTE-elements. 

3. Hydrodynamical Models 

3.1. EXPLOSIONS OF MASSIVE WHITE DWARFS 

A first group consists of massive carbon-oxygen white dwarfs (WDs) with a 
mass close to the Chandrasekhar mass which accrete mass through Roche- 
lobe overflow from an evolved companion star (Nomoto &: Sugimoto 1977; 
Nomoto 1982). The explosion is triggered by compressional heating. The 
key question is how the flame propagates through the white dwarf. Several 
models of SNela have been proposed in the past, including detonation (Ar- 
nett 1969; Hansen & Wheeler 1969), deflagration (e.g. Nomoto, Sugimoto 
&; Neo 1976) and the delayed detonation model, which assumes that the 
flame starts as a deflagration and turns into a detonation later on (Khokhlov 
1991, Wooslcy k Weaver 1995, Yamaoka ct al. 1992). 

Our sample includes detonations (DETl/2), deflagrations (W7, DFl, 
DFlmix), delayed detonations (N21/32, M35-39; DD13-27) and pulsating 
delayed detonations (FDD 1-9). The deflagration speed is parameterized as 
Ddef = 1 where is the local sound velocity ahead of the flame and a is 
a free parameter. The speed of the detonation wave is given by the sound- 
speed behind the front. For delayed detonation models, the transition to a 
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TABLE 1. Some quantities (see text) are given for detonations (DET1/DET2, asterix), de- 
flagrations (W7, DFl, DFlmix, circles), delayed detonations (N21, N32, M35-M312, DDlSc 
- 27c, bullets), pulsating delayed detonations (PDDla-9, black triangles), envelope models 
(DET2env2/4, open triangles), and helium detonations (HeD2-12, open asterix). The C/0 
ratio has been assumed to be 1:1 unless otherwise quoted in brackets after the name. The 
initial metallicity for > 20 is assumed to be solar, but for DD24c, DD25c, DD26c and 
DD27c for which 1/3, 3, 0.1 and 10 times solar abundances arc used, respectively. Otherwise, 
these models are identical to DD21c. For the hcHum detonations and envelope models, the 
mass is of the C/O core, and of the He-layers or of the extended CO-envelope, respectively, 
are given separately. In HeD7, no central C-ignition is triggered. 
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Figure 4- Abundances as a function of the expansion velocity for three delayed detona- 
tion models (see table 1). Since the burning time scales to NSE or partial NSE are shorter 
than the hydrodynamical time scales for all but the very outer layers, the final products 
depend mainly on the density at which burning taJfes place. With decreasing transition 
density, lesser ^^Ni is produced and the intermediate mass elements expand at lower ve- 
locities because the later transition to a detonation allows for a longer pro-expansion of 
the outer layers (DD21c vs. DD13c). Siniilarily, with increasing C/0 ratio in the progen- 
itor, the specific energy release during the nuclear burning is reduced (DD23c vs. DD21c) 
and the transition density at the burning front reached later in time, resulting in a larger 
preexpansion of the outer layers. This may allow to determine the main sequence mass 
of the progenitor. 



detonation is given by another free parameter ptr- When the density ahead 
of the deflagration front reaches pu- , the transition to a detonation is forced 
by increasing a to 0.5 over 5 time steps bringing the speed weh above the 
Chapman-Jouguet theshold for steady deflagration. For pulsating delayed 
detonation models, the initial phase of burning fails to release sufficient 
energy to disrupt the WD. During the subsequent contraction phase, com- 
pression of the mixed layer of products of burning and C/0 formed at the 
dead deflagration front would give rise to a detonation via compression and 
spontaneous ignition (Khokhlov 1991). In this scenario, ptr represents the 
density at which the detonation is initiated after the burning front dies out. 
Besides the description of the burning front, the central density of the WD 
at the time of the explosion is another free parameter. For white dwarfs 
close to the Chandrasekhar limit, it depends sensitively on the chemistry 
and the accretion rate M at the time of the explosion. In our recent study 
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Figure 5. Abundances for three delayed detonation models (DD21c, DD24c & DD27c) 
with identical parameters but different initial metallicity (Z > 20) corresponding to 1, 
1/3 and 10 times solar. The overall density, velocity and chemical structure remains 
uneffected but the isotopic composition changes drastically. In particular, the final Fe 
abundance differ by more than two orders of magnitude in the outer layers and may be as 
high as 7 % which clearly shows up in synthetic spectra. Whether this provides a natural 
explanation of the spectral peculiarities of SN1991T is under investigation. 



on SN1994D, evidence was found that the models with a somewhat smaller 
central density provide better agreement with both the observed spectra 
and LCs. Therefore, the grid of PDD models has been extended. 



3.2. MERGING WHITE DWARFS 



The second group of progenitor models consists of two low-mass white 
dwarfs in a close orbit which decays due to the emission of gravitational 
radiation and this, eventually, leads to the merging of the two WDs (e.g. 
Iben &; Tutukov 1984). After the initial merging process, one low density 
WD is surrounded by an extended envelope (Hachisu, Eriguchi &; Nomoto 
1986ab, Benz, Thielemann & Hills 1989). This scenario is mimicked by our 
envelope models DET2env2...6 in which we consider the detonation in a 
low mass WD surrounded by a compact envelope between 0.2 and 0.6 Mq. 
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Figure 6. Monochromatic V LCs of normal bright supernovae of deflagration (W7), 
classical delayed detonation (M36), Helium detonation (HeDlO), envelope (DET2env2), 
pulsating delayed detonation (PDD3) models, and of subluminous SNela classical delayed 
detonation (M39), Helium detonation (HeD6), pulsating delayed detonations(PDD5). 



3.3. EXPLOSIONS OF SUB-CHANDRASEKHAR MASS WHITE DWARFS 

Another class of models - double detonation of a C/ O- WD triggered by det- 
onation of helium layer in low-mass WDs - was explored by Nomoto (1980), 
Woosley Sz Weaver (1980), and most recently by Woosley and Weaver (1994, 
hereafter WW94) and Hoflich & Khokhlov (1995). This scenario was also 
suggested for the explanation of subluminous Type la (WW94). Note that 
the explosion of a low mass WD was also suggested by Ruiz-Lapuente et 
al. (1993) to explain subluminous SNe la but the mechanism for triggering 
of the central carbon detonation was not considered. 

For ease of comparison, wc have used parameters close to those sug- 
gested in WW94. To prevent repetition, we refer to the latter work for 
a detailed discussion of this class of models. Helium detonations show a 
qualitatively different structure in comparison to all models with a Chan- 
drasekhar mass WD . The intermediate mass elements are sandwiched by 
Ni and He/Ni rich layers at the inner and outer regions, respectively. Gen- 
erally, the density smoothly decreases with mass because partial burning 
produces almost the same amount of kinetic energy as the total burning, 
but a moderate shell-like structure is formed just below the former Helium 
layers. Observationally, a distinguishing feature of this scenario is the pres- 
ence of Helium and Ni with expansion velocities above 11,000 to 14,000 
km s~^. Typically 0.07 to 0.13 Mq of Ni are produced in the outer layers, 
mainly depending on the mass of the Helium shell. Recently, Benz (1995, 
private communication) suggested that the explosion of the C/0 core may 
be triggered directly by the in-going shock front if Helium is ignited some- 
what above the He/C-0 interface. This, however, should hardly change the 
chemical structure because the flame propagates as a detonation. 
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Figure 7. Ratio Q between the bolometric luminosity and the energy release by gamma 

decay at ttoi as a function of the ''^Ni mass (left) and rise time to bolometric maximum 
tboi- The different symbols correspond to different explosion scenarios (see Table 1). 
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Figure 8. Absolute visual brightness V as a function of the ^^Ni mass (left) and intrinsic 
color B-V as a function of V(right). 



4. Comparison with Observations 

The different explosion scenarios can generally be distinguished based on 
differences in the slopes of the early monochromatic LCs (Fig. 6) and the 
expansion velocities indicated by the spectra (Fig. 9 & Hofiich, this volume). 
For all models with a ^^Ni production > QAMq, My ranges from -19.1 to 
— 19.7™. As a general tendency, the post-maximum declines are related to 
My, but there is a significant spread in the decline rate among models with 
similar brightess. For all models but the Helium detonations, the colors 
become very red for small Mjvi (Fig. 8). 

As a general trend, the maxima of subluminous supernovae are more 
pronounced. The reason becomes obvious if we consider the following re- 
lated point. In the literature (Arnett et al. 1985), it is often assumed that 
the bolometric luminosity L^q; is equal to the instantaneous release of en- 
ergy by radioactive decay L-y, i.e. Q = L^oi/E^ « 1, independent from the 
model. We find that Q depends on the model because the opacities vary 
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strongly with temperature. Typically, Q ranges from 0.7 to 1.4 (Fig. 7). 
For models with fast rising LCs, the opacity stays high and the photo- 
sphere recedes mainly by the geometrical dilution of matter. For models 
with a slower rise time or little ^^Ni, the opacity drops strongly at about 
maximum light. The photosphere receds quickly in mass, and thermal en- 
ergy can be released from a larger region. Because no additional energy 
is gained, the energy reservoir is exhausted faster, and the post-maximum 
decline becomes steeper. 

The distinguishing features of Helium detonations are blue maxima both 
for normal and subluminous supernovae, a rapid increase of the luminosity 
caused by ^^Ni heating, followed by a phase of slow rise to maximum light 
and a fast post-maximum decline. The fast decline is caused by the rapid 
increase of the escape probability of 7 photons originating in the outer ^^Ni 
and the short life time of ^^Ni. The early synthetic spectra are dominated 
by high velocity Ni and the absence of Si which is underabundant by more 
than two orders of magnitude compared to delayed detonation models. 
Note that a significant amount of the line emission at late phases should be 
powered by ^^Co at high velocities. In the subluminous models, about 1/2 
to 1/4 of the radioactive material is ejected at high velocities. Accordingly, 
late-time spectra should show broad Co and Fe lines. 

5. Comparison between Observation and Models Predictions 

We use a quantitative method for fitting data to models based on Wiener 
filtering (Rybicki and Press, 1995). The reconstruction technique is applied 
to the standard deviation from the theoretical LC to avoid problems with 
measurements distributed unevenly in time. By minimizing the error, the 
time of the explosion, the distance, and the reddening correction can be 
determined. For details, see Hoflich & Khokhlov (1995). 

Observed monochromatic LCs and spectra of 27 SNe la are compared 
with theoretical models (Table 2, Fig. 9 & Hoflich, this volume). According 
to our results, normal bright, fast SNela (e.g. SN 1971G, SN 1994D) with 
rise times up to 15 days (17 days) for the blue (visual) LC can be explained 
by delayed detonation with different densities ptr for the transition from 
a deflagration to a detonation. For PDDs, the density ptr stands for the 
density at which the detonation starts after the first pulsation. Typically, 
Ptr is about 2.5 lO''' g cm~^. Central densities of the initial WDs range 
from 2.1 to 3.5 10^ g cm~^. As a tendency, models at the lower end of 
this range tend to give somewhat better fits. This may be explained by a 
high accretion rate, some variation in the chemical composition or by an 
additional trigger mechanism for the explosion. We note that the classical 
deflagration W7 (Nomoto et al. 1984) provides similar good fits in several 
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TABLE 2. List of observed SNela used in our sample. Columns 2 to 7 give the parent galaxy, 
its type, peculiar velocity Vz according to the MERCG catalogue (Kogoshvili 1986), the 
distance modulus of the host galaxy and the color excess Eb-v of the SNe la as determined 
from our models. In column 8, we give names of those theoretical models (see Table 1) which 
can reproduce the observed LCs. Models in bracketts do not fuUfil our criterion but are 
close. These models have been excluded for M-m. Models marked with ^ can be ruled 
out because the observations indicate Si at high velocities incompatible with the chemical 
profiles. Note that the new delayed detonation models have not been included in this list. 
We found some evidence for a relation between the type of the explosion and of the host 
galaxy. In elliptical galaxies, SN la with shell-like structures are highly favored. These may 
be understood within the pulsating and merging scenarios. Because ellipticals consist on 
low mass stars only, this may provide a hint to the progenitor evolution; however, a larger 
sample of observations is needed to confirm this trend. 
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cases because its structure resembles those of DD models but it has some 
problems with the high velocities of Si lines in SN1994D (Hoflich 1995). 




Figure 9. Monochromatic LCs of SN 1989B band compared with the calculated LC of 
the delayed detonation model M37. The 2 cr-error ranges are given by Wells et al. (1994). 
In the lower left plot, the photosphcric velocity of the delayed detonation model M37 is 
given in comparison with the expansion velocities implied by the line shifts of Mg II and 
Sill. The dotted line is the velocity at which, in M37, the Si abundance has dropped to 
10 % of its maximum value. The 'leveling off' of the line shift is produced by the chemical 
profile and not due to a change of the density profile. Note that detailed spectral synthesis 
is required to archieve a better accuracy or to test the ionization balance (Hoflich 1995) . 



The "standard" explosion models are unable, however, to reproduce rise 
times to blue (visual) maximum longer than 15 days (17 days), provided 
the progenitor is a C/0 WD of about 1.2 to IAMq. In fact, slow rising and 
declining LCs have been observed (e.g. SN 1990N, SN 1991T) which require 
models with an envelope of typically 0.2 to QAMq. The envelope can be 
produced during a strong pulsation or during the merging of two WD. The 
lower value should not be regarded as a physical limit, because it is likely 
that a continuous transition from models with and without an envelope 
exists. Note that a unique feature of models with massive envelopes are 
very high photospheric expansion velocities {vph ~ 16,000 km/s) shortly 
before maximum light, which drop rapidly to an almost constant value 
between 9000 and 12, 000 km/s. This "plateau" in Vph lasts for 1 to 2 weeks 
depending on the envelope mass (Khokhlov et al. 1993). In fact, there is 
some evidence for the plateau in Vph from the Doppler shift of lines of SNe la 
with a slow pre-maximum rise and post-maximum decline (e.g. SN 1984A, 
SN 1990N, Miiller & Hoflich 1994). 

Strongly subluminous supernovae (SN1991bg, SN1992K, SN1992bc) can 
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be explained within the framework of pulsating delayed detonation mod- 
els with a low transition density. In particular, the models become sys- 
tematically redder and the post-maximum decline becomes steeper with 
decreasing brightness in agreement with observations. However, we do not 
get unique relations between these different quantities. The evolution of the 
photospheric expansion velocity Vph (Hoflich et al. 1995) and, in particu- 
lar, its steady decline, is consistent with observations. We must also note 
that we need a rather high intrinsic reddening for SN1992K and SN1991bg. 
Whether this can be explained by selective line blanketing, dust formation, 
or foreground clouds, or a combined effect is under investigation. The latter 
possibility must be regarded as unlikely, because, at later phases, the colors 
of these two supernovae are close to those of bright SNe la. Note that some 
of the 'classical' delayed detonation models (M39) also produce strongly 
subluminous LCs, but these do not fit any of the measurements. 

Our Helium detonation models arc rather Tinsuccessful in reproducing 
observations, mainly due to the rather steep post-maximum declines for 
normal bright supernovae, and Si lines observed at higher expansion ve- 
locities than compatible with the chemical profiles. At early times, our 
synthetic spectra are dominated by Fe-group elements and lack strong 
Si features because of the low abundances of Si (< 10~^" ~^). For sub- 
luminous supernovae, the blue color at maximum light and the strong IR- 
maximum are both in contradiction to the observations. Quantitatively, 
multi-dimensional effects may alter the LCs mainly due to an increased 
escape probability for photons. However, the basic features of the LCs and 
spectra are not expected to change because they are inherent to the outer 
^^Ni. The energy required to push the entire ^^Ni to sufficient high veloc- 
ities (< 16, OOOA;m/sec) is well in excess of the entire energy production. 

Our findings with respect to the explosion scenario can be concluded as 
follows. Models with masses close to the Chandrasekhar limit provide the 
best agreement with the observations. Delayed-detonation and deflagration 
models similar to W7 and pulsation or merging scenarios are required. 
Based on the sparse observations, we cannot exclude the mechanism of 
Helium detonations. In particular, more early time spectra are crucial. 

Pinto's group end up with a different result. Based on their calculations 
of the delayed detonation model DD4, they exclude this scenario. Their B 
and V LCs show rise times in excess of 21 days whereas both the U and 
bolometric LCs peak around day 15 and, at maximum light, their optical 
spectra show Co lines in strong emission (talk by R. Eastman). However, 
according to Pinto (priv. communication), their Helium detonations pro- 
vide both excellent fits to spectra and LCs without encountering our Si 
problem. We find Hg to be 67 ± 9km/secMpc within a 95 % confidence 
level. This value agrees well with our previous analysis based on a subset of 
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Figure 10. Visual LC of SN 1988U in comparison with the calculated LC of the delayed 
detonation model M36. Both the V and R colors are used. The dashed line marks the 
time when the photospheric velocity of the models corresponds to the observed line shift 
(thin line). The right and left plots correspond to different time shifts and both allow for 
a reasonable reproduction of the LC. This would correspond to an uncertainty of 0.5™ 
in the brightness. However, the right plot can be ruled out by the spectrum. This clearly 
demonstrates the importance of simultaneous analysis of both LCs and spectra. 



observations and models (66 it IQkm Mpc~^ sec"^ , Miiller &; Hofiich 1994). 
A strong variation of the local value can be ruled out at least on scales 
below 200 Mpc. From SN1988U, the deceleration parameter qo is 0.7 it 1.0 
(2(t). Better limits can we expected in the near future when more and more 
distant SNela become available (Perlmutter et al. 1995). 



6. Distance determinations, Hg and qo 

Based on our LCs, we have also determined the individual distances of the 
parent galaxies of the analyzed SNe la (Table 2). Our method does not 
rely on secondary distance indicators and allows for a consistent treatment 
of interstellar reddening and the interstellar rcdshift. The advantages of a 
consistent inclusion of information from the spectra becomes striking for 
SN1988U (Fig. 10). 

Other determinations of Hq which are based on independent, purely 
statistical methods and primary distance indicators. It may be encouraging 
that the result of different SN la based methods agree if SN la are not 
treated as as standard candles. Hamuy et al. (1995) found 65ib5, Riess et al. 
(1995) give 67±5, and Fisher et al. (1995), Nugent et al. (1995b) and Branch 
et al, (1995) get a values of 60 ± 10, 60 ± 12 and 58 it 7, respectively. Prom 
our models, both the empirical relations between My/(iM(15) like-relations 
and the ansatz to deselect subluminous SNela seems to be justified, but we 
expect an individual dispersion of a; 20% (Hamuy et al. 1995b). 
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Figure 11. Hubble values H are shown based on individual distances (see table 2). 
SN1988U at v=91500 km/sec gives = 64 ± 10 [km/sec Mpc]. 
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